home *** CD-ROM | disk | FTP | other *** search
/ Shareware Grab Bag / Shareware Grab Bag.iso / 005 / miscstat.arc / EXPREG.BAS < prev    next >
Encoding:
BASIC Source File  |  1986-04-25  |  1.2 KB  |  49 lines
  1. 1 'EXPonential REGression
  2. 2 KEY OFF:SCREEN 0,1,0:WIDTH 40:COLOR 15,1,9:CLS
  3. 10  PRINT "EXPONENTIAL REGRESSION"
  4. 20  PRINT 
  5. 30  PRINT "NUMBER OF KNOWN POINTS";
  6. 40  INPUT N:IF N=0 THEN 410
  7. 42 IF N<1 THEN 30
  8. 50 J = 0
  9. 60 K = 0
  10. 70 L = 0
  11. 80 M = 0
  12. 90 R2 = 0
  13. 100  FOR I = 1 TO N
  14. 110  PRINT "X,Y OF POINT";I;
  15. 120  INPUT X,Y
  16. 130 Y =  LOG (Y)
  17. 140 J = J + X
  18. 150 K = K + Y
  19. 160 L = L + X ^ 2
  20. 170 M = M + Y ^ 2
  21. 180 R2 = R2 + X * Y
  22. 190  NEXT I
  23. 200 B = (N * R2 - K * J) / (N * L - J ^ 2)
  24. 210 A = (K - B * J) / N
  25. 220  PRINT 
  26. 230  PRINT "A =" EXP (A)
  27. 240  PRINT "B ="B
  28. 250 J = B * (R2 - J * K / N)
  29. 260 M = M - K ^ 2 / N
  30. 270 K = M - J
  31. 280  PRINT 
  32. 290 R2 = J / M
  33. 300  PRINT "COEFFICIENT OF DETERMINATION (R^2): ";R2
  34. 310  PRINT "COEFFICIENT OF CORRELATION: "; SQR (R2)
  35. 320  PRINT "STANDARD ERROR OF ESTIMATE: "; SQR ( ABS (K / (N - 2)))
  36. 330  PRINT 
  37. 340  PRINT "INTERPOLATION: (ENTER 0 TO END PROGRAM)"
  38. 350  INPUT "X = ",X
  39. 370  IF X = 0 THEN 410
  40. 380  PRINT "Y = "; EXP (A) *  EXP (B * X)
  41. 390  PRINT 
  42. 400  GOTO 350
  43. 410 PRINT:LOCATE 24,1:PRINT"PRESS <SPACE> TO RUN AGAIN";
  44. 420 LOCATE 25,1:PRINT"PRESS <ESC> TO EXIT";
  45. 430 K$=INKEY$:IF K$="" THEN 430
  46. 440 IF K$=" " THEN RUN
  47. 450 IF K$<>CHR$(27) THEN 430
  48. 10000 RUN"EXPREG
  49.